Magnetically- induced ferroelectricity in orthorhombic manganites: 
microscopic origin and chemical trends 
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The microscopic origin of the magnetically-driven ferroelectricity in collinear AFM-E orthorhom- 
bic manganites is explained by means of first-principles Wannier functions. We show that the 
polarization is mainly determined by the asymmetric electron hopping of orbitally-polarized Cg 
states, implicit in the peculiar in-plane zig-zag spin arrangement in the AFM-E configuration. In 
ortho-HoMnOs, Wannier function centers are largely displaced with respect to corresponding ionic 
positions, implying that the final polarization is strongly affected by a purely electronic contribution, 
at variance with standard ferroelectrics where the ionic-displacement is dominant. However, the fi- 
nal value of the polarization is the result of competing effects, as shown by the opposite signs of the 
contributions to the polarization coming from the Mn Cg and t2g states. Furthermore, a systematic 
analysis of the link between ferroelectricity and the spin, orbital and lattice degrees of freedom in 
the manganite series has been carried out, in the aim of ascertaining chemical trends as a function of 
the rare-earth ion. Our results show that the Mn-O-Mn angle is the key quantity in determining the 
exchange coupling: upon decreasing the Mn-O-Mn angle, the first- (second-) nearest neighbor ferro- 
magnetic (antiferromagnetic) interaction decreases (remains constant), in turn stabilizing either the 
AFM-A or the AFM-E spin configuration for weakly or strongly distorted manganites, respectively. 
The Mn eg contribution to the polarization dramatically increases with the Mn-O-Mn angle and 
decreases with the "long" Mn-0 bond length, whereas the Mn t2g contribution decreases with the 
"short" Mn-0 bond-length, partially cancelling the former term. 

PACS numbers: 75.47.Lx,75.80.+q, 75.50.Ee, 77.80.-e 



I. INTRODUCTION 

Multiferroics are attractive multifunctional materials 
where magnetism and ferroelectricity coexist and are gen- 
erally coupled. In particular, orthorhombic rare-earth 
manganites i^MnOa represent an important class of "im- 
proper ferroelectrics" [l| where electric dipoles are in- 
duced by a frustrated magnetic ordering. Within this 
family of compounds, TbMnOa and DyMnOa, in their 
noncollinear magnetic phases, have been experimentally 
shown to behave as multiferroics showing a weak polar- 
ization (P < 0.1 /iC/cm^) and a spin- flop ferroelectric 
transition. [2] Moreover, it has been recently predicted 
that relatively strong ferroelectricity occurs in the E- 
type antiferromagnetic (AFM) phase of i^MnOa through 
a model study where the double-exchange interaction be- 
tween Mn-d orbitals is proposed as a driving force for 
polar atomic displacements. [3] In that case, the electric 
polarization P is not related to Dzyaloshinskii-Moriya 
interaction, which is expected to generate much lower 
P. [4, 5, 6] Following the model study, our previously 
reported ab-initio calculations [7] have confirmed that 
AFM-E HoMnOa indeed shows a high ferroelectric po- 
larization (P ~ QjiC /cm?)] it was shown there that the 
AFM-E spin ordering, which breaks the space-inversion 
symmetry, is such that ferroelectric dipoles arise even 
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without atomic displacements. 0] From the experimental 
point of view, ferroelectricity was detected in policrys- 
talline ortho-HoMn03,[8] the magnitude of polarization 
being however much smaller than what is theoretically 
predicted and showing a strong dependence on the mag- 
netic field below the ordering temperature of the Ho, sug- 
gesting their involvement in the development of P. The 
reason behind the disagreement between theory 0, 0] and 
experiments [8] is still under debate. 

The microscopic origin of the multiferroism in ortho- 
i^MnOa is tightly linked to the lattice degree of freedom: 
a small ionic radius of the R atom holds a key role in af- 
fecting i) structural properties, such as both the JT dis- 
tortion in the Mn02 plane and the GdFeOs-like tilting of 
MnOe octahedron, ii) spin configuration determined by 
Mn-d super-exchange interaction Jij , Hi) orbital ordering 
stabilized by Jahn- Teller (JT) distortion, and iv) hopping 
integral in terms of double exchange interaction. Since 
these aspects are all mutually combined, first we will dis- 
cuss the case of AFM-E HoMnOa, and we will explain, 
from a microscopic quantum-mechanical point of view, 
the origin of the magnetically-induced ferroelectric po- 
larization in terms of the spin, orbital and lattice degrees 
of freedom. The analysis based on Wannier functions 
represents a novel development with respect to our pre- 
vious study on HoMnOs.Q Moreover, in the second part 
of the paper, we will discuss about structure, magnetism, 
orbital ordering, hopping integrals and ferroelectricity - 
as well as the links between them - along the manganites 
series, with the main aim of identifying chemical trends 
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as a function of the R ion. 



II. STRUCTURAL AND COMPUTATIONAL 
DETAILS 

The unit ceh in orthorhombic i^MnOa shows strong 
distortions with respect to the ideal cubic perovskite. 
Whereas the Pbnm setting is used in some other refer- 
ences, the standard Pnma orientation is adopted in this 
paper, i.e. we choose h as the longest axis. Although 
the hexagonal non-perovskite phase is more stable than 
the orhorhombic phase in i^MnOa for R smaller than 
Tb, the transition to the metastable orthorhombic phase 
can be obtained by high-pressure synthesis. [ol. fiol. Illl. [l2j 
We have considered five kinds of AFM configurations: 
AFM-A, C, G in 20 atoms/unit ceh and AFM-E and 
AFM-E* in 40 atoms/unit cell. We recall that, according 
to the standard Wollan-Koehler notation[13], the AFM- 
A shows FM (AFM) intraplanar (interplanar) coupling; 
the AFM-C shows FM (AFM) interplanar (intraplanar) 
coupling; the AFM-G shows AFM in both intra- and in- 
terplanar coupling. The AFM-E shows in-plane FM zig- 
zag chains antiferromagnetically coupled to the neighbor- 
ing chains; the interplanar coupling is also AFM (for the 
AFM-E spin arrangements in the Mn02 plane, we refer to 
the following Fig. [2]). We denote by AFM-E* the spin- 
configuration showing the same in-plane spin arrange- 
ment as AFM-E, but with an interplanar FM coupling. 
Note that the AFM-A spin arrangement shows the space- 
inversion as symmetry operation, at variance with the 
AFM-E spin-configuration which shows non centrosym- 
metric Pmn2i symmetry. 

As far as the electronic structure is concerned (espe- 
cially for the density of states projected on the Mn atom 
as well as orbital-ordering, see below), it is useful to de- 
fine - in addition to the "global" X, F, Z orthorhombic 
frame - a "local" frame, specific to each J ahn- Teller- like 
distorted MnOe octahedron, obtained by choosing x, z 
along the middle, short and long Mn-0 axis, respectively. 
In this local frame, the orbital ordered Mn-eJ state, which 
is often expressed as (3x^ ~ ^^)/ (3^^ — ^'^)^ is described as 
3z^ — r^ (denoted as z'^ hereafter for simplicity). The sim- 
ulations were performed by using density-functional the- 
ory (DFT) and the Perdew-Becke-Erzenhof (PBE) ver- 
sion of the generalized gradient approximation (GGA) to 
the exchange-correlation potential. [M] The calculations 
were done with two program codes according to different 
purposes. 

The expensive calculations for the structural op- 
timization of the atomic structure and Berry phase 
for the AFM-A and AFM-E phases were done with 
"Vienna Ab initio Simulation Package (VASP) pro- 
gram code" ,[l5| where projector- augmented- wave (PAW) 
pseudo-potential is used. The plane wave cutoff is set 
to 500 eV in energy. We used 8 special k-points (di- 
vided as 4 X 3 X 4) in 1/8 irreducible Brillouin zone (IBZ) 
for the A- AFM phase and 4 special k-points (divided as 



4 




-7 - 



Y r X s r z R 

FIG. 1: Energy bandstructure of optimized AFM-E HoMnOs- 
Due to the AFM-E state, only the up-spin channel is plotted. 

denotes the Fermi energy and it is set as zero of the energy 
scale. 



2 X 3 X 4) in 1/4 IBZ for the E-AFM phase according to 
the Monkhorst-Pack scheme. The Berry phase was cal- 
culated by integrating over six k-point strings parallel to 
the c axis, each string containing 6 k-points. 

Calculations requiring a higher precision {i.e. as 
for the total energy differences for the stability of 
different spin-configuration) and the construction of 
Maximally-localized Wannier functions (WFs) [l^ were 
done with the FLEUR codefT?] which is based on the 
full-potential linearized augmented planewave (FLAPW) 
formalism. [18] Muffin-tin radii were set to 2.5, 2.0 and 
1.5 a.u. for i?, Mn and O atoms, respectively, where the 
wavefunction cutoff was chosen as 3.8 a.u.~^. The po- 
tential was converged with 24 special k-points and den- 
sity of states was calculated with 192 k-points within the 
tetrahedron method. For AFM-E and AFM-E* phase, 
12 k-points were used according to the doubled unit cell. 
The Wannier function calculation, whose procedure was 
recently implemented in the FLEUR code [l^ , was done 
with 512 k-points (divided as 8 x 8 x 8). R-bs and bp 
states were treated as local orbitals. 

The localized i?-4/ electrons were assumed as core elec- 
trons: "frozen core" within the VASP code and "open 
core" within the FLEUR code, [20] where the spin mo- 
ment is maximized due to Hund's rule. Irrespective of 
these approaches, the 4/ states lie deep in energy (a few 
eV below the Fermi energy) and they are almost com- 
pletely undispersed, so that they do not affect other va- 
lence states. 



III. MICROSCOPIC ORIGIN OF THE 
FERROELECTRIC POLARIZATION IN HoMnOs 

Let us start our discussion by focusing on the case 
of AFM-E HoMnOa, which was previously found by 
some of us to show a large polarization along the in- 
plane c-axis due to non-centrosymmetric collinear spin- 
arrangement. [7] In this work, we perform a careful anal- 
ysis in terms of Wannier functions which explains the 
microscopic mechanism at the basis of the final ferroelec- 
tricity. The structural parameters used in our calcula- 
tions are summarized in a later section, together with 
the parameters in other i^MnOa compounds (cfr Table 
IIV A II below). The magnetically-induced spontaneous 
polarization was calculated by using i) the point charge 
model (PCM) where each ion has been given its nominal 
charge (Ho:+3, Mn:3+, 0:2-); ii) the Berry phase (BP) 
method implemented in the VASP code and in) the Wan- 
nier function (WF) method implemented in the FLEUR 
code.[2lj We recall that, in the first approach only the po- 
sitions of the anions and cations are considered, whereas 
in the two latter (2^, [23] quantum-mechanical treatments, 
the self-consistent electronic structure is fully taken into 
account. 

To calculate the polarization within BP, an adia- 
batic path from AFM-A to AFM-E phase is assumed, 
in such a way that the direction of Mn spins are pro- 
gressively rotated from an in-plane FM to a zig-zag- 
like arrangement. [22] The ionic contribution from core 
electrons and protons are added to BP calculated for 
the fully occupied valence states. In Fig|2] and Table 
H] we show the displacements of the atoms in the non- 
centrosymmetric AFM-E spin-configuration with respect 
to the centrosymmetric AFM-A structure. These quan- 
tities will be needed in the discussion reported below. 

By using WFs, one can decompose the total polariza- 
tion into contributions coming from each set of orbit als. 
For clarity, we show in FigHlthe AFM-E HoMnOs insu- 
lating bandstructure where the relevant states are high- 
lighted: Mn-e^, Mn-t2g and 0-p orbitals. We have pro- 
jected these three groups of occupied eigenstates into 
real-space basis separately and "maximally-localized" 
them to obtain the corresponding WFs. The contribu- 
tion from deeper occupied valence states (such as 0-25 
and i?-5s, -5p states) is neglected in the WF approach. 
Here, the total polarization is the sum of the displace- 
ments of the centers of each WF from the position of the 
corresponding ion plus PCM contribution. For the de- 
tails of the construction and interpretation of WFs, see 
RefHi. 

According to the different approaches, our estimated 
values for the polarization in AFM-E HoMnOa are Pbp 
= -6.14 /iC/cm^, PwF = -5.61 /iC/cm^ and PpcM = - 
2.09 /iC/cm^ along Z axis. The values of P along X 
and Y axis are both negligible in each approach. The 
large difference of the values of Pwf (Pbp) from PpcM 
is remarkable. 




FIG. 2: Atomic displacements in FE HoMnOs, as obtained 
by the difference of atomic coordinates in optimized AFM-E 
and optimized AFM-A spin configuration (length of arrows in 
arbitrary units) . A schematic representation of the Cg orbital- 
ordering and the direction of the ferroelectric polarization are 
also shown. 



TABLE L Atomic displacements from centrosymmetric AFM- 
A phase in ferroelectric AFM-E HoMnOs (A). 0(1)^ and 
0(2)P (0(3)^P and 0(4)^^) are connected by C2Z symmetry 
operation. The used notation of atoms is shown in Fig. 2. 





AX 


Ay 


AZ 


\dr\ 


Mn(0) 
0(1)P 
0(3)"P 


0.036 
-0.006 
0.020 


0.003 
0.009 
0.010 


-0.010 
0.016 
0.009 


0.037 
0.020 
0.024 



A. Wannier Representation of Anomalous 
Contributions 

The difference between polarization Pwf (^bp) and 
PpcM is referred to as the anomalous contribution[2E\ 
and it is analyzed by means of Wannier functions. Table 
im shows the orbital-decomposed polarization obtained 
by the displacement of the center of WF (WFC) with 
respect to ionic positions for each set of bands. By defi- 
nition, the sum of the WFC displacements gives the dif- 
ference between Pwf and PpcM, which can be regarded 
as the difference between the dynamical charge Z* and 
static (nominal) charge Z. In other words, it quantifies 
the effects of going from a more ionic charge distribution 
to more covalent bonding state. By means of the WF 
analysis, not only one can obtain relevant information 
on which orbital causes the anomalous contribution, but 
also a real-space picture of the polar orbital states. 

Let us first focus on the e^-like d^^ and dyz orbitals 
which, together with 0-p orbitals, mostly contribute to 
polarization (see main displacements along the polar c 
axis inTable. 

First, consider the e^-orbital at Mn(0)-up site sur- 
rounded by four O ions and four second-neighbor Mn 
ions in the ac plane (Figl3]a). Because of the JT-derived 
orbital ordering, the Cg-.d^^ orbital of Mn(0) ion is point- 
ing towards the two O ions along long axis {z in the 
"local" frame). Therefore, it makes a a bonding with 
0(1)P- and 0(3)^P-p2 orbitals. In addition, when con- 
sidering the AFM-E spin configuration, the up-spin elec- 
tron of the e^-orbital of the Mn(0) ion can hop only onto 
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(a) Mn-Cg : 




FIG. 3: Isosurface of WFs for Mn-d states centered at Mn(0) 
site in AFM-E HoMnOs: a) eg states and b) t2g states. Green 
arrows indicate the displacement of the centers of WFs from 
the Mn atomic position. Here, the superscript p (ap) of O 
denotes oxygen ion between parallel (antiparallel) spin of Mn 
ions. 




FIG. 4: Isosurface of WFs of 0-2p state. Only four typical 
WFs in the ac plane are shown, a) -py, b) O'^-pz, c)0^^-py^ 
d)0^^-pz (up spin state). Arrows indicate the displacement 
of the centers of WFs from the O atomic position. 



TABLE II: WFC displacements from atomic coordinates in 
AFM-E HoMnOa (A). The orbitals are denoted in a local 
frame (ximiddle, yishort, zilong axis for Mn, x:antibonding 
inter-plane, y:antibonding in-plane, z:bonding for O). Only 
up spin contribution is shown. 









AX 


Ay 


AZ 


\dr\ 


Mn(0) 


eg : 


z' 


-0.169 


0.044 


0.188 


0.257 




t2g 


: xy 


0.021 


-0.008 


-0.017 


0.028 




t2g 


: yz 


-0.049 


0.018 


-0.146 


0.155 




t2g 


: zx 


0.013 


0.007 


-0.008 


0.016 


0(1)P 


Px 




-0.158 


0.071 


0.200 


0.265 




Py 




-0.110 


0.048 


0.209 


0.241 




Pz 




0.004 


0.037 


0.108 


0.115 


0(3)"P 


Px 




0.020 


-0.001 


-0.026 


0.033 




Py 




-0.025 


-0.013 


-0.038 


0.047 




Pz 




-0.171 


-0.016 


0.048 


0.178 




Mn(l)f 0(1) Mn(0)f 0(3) Mn(3)| 



FIG. 5: Schematic representation of the mechanism which 
causes the microscopic polarization of Mn-e^ orbital due to 
the asymmetric electronic hopping. Bottom sketch: Hop- 
ping is allowed only from Mn(0) to Mn(l) and not to Mn(3): 
this determines the direction of WFC (empty arrow) . To in- 
crease hopping, Mn(0) moves "right" (see filled arrow), giving 
a lower weight of WF on Mn(0) but larger on Mn(l) [repre- 
sented as dashed vs solid (before vs after displacement) Gaus- 
sians on the atoms] . Also shown is the resulting P (see arrow) . 
Top sketch: Resulting atomic displacement (shown as filled 
arrows) aimed at increasing 0^. 



Mn(l) (up spin) site via the 0(1)^ atom but not on the 
Mn(3) (down spin) site. This asymmetric hopping gives 
the considerable deviation of the WFC from the atom 
position (0.26 A) whose direction is shown as a green ar- 
row in FiglSti. Recall that the hopping integral between 
Mn-e^ orbitals depends on both the Mn-0 bond length 
d and the Mn-O-Mn angle (/), the hopping energy being 
therefore optimal for larger (j). As a consequence, in or- 
der to increase the hopping, the Mn(0) ion is expected 
to move toward 0(3)^^ to increase (j)^ between parallel 
Mn spins, i.e. along a direction which is opposite to 
the WFC displacement. In the same aim of increasing 
hopping, 0(1)P moves perpendicularly to the Mn^-Mn^ 
bonding. Similarly, the 0(3)^^ atom is displaced perpen- 
dicularly to the Mn^-Mn^ bonding (see rough schema- 
tization of atomic displacements shown in Fig|5]). This 
is consistent with what is shown in Figl2] and Table [D 
The ions therefore move to enhance the polarization in- 
duced by the asymmetric hopping of Cg -orbital electrons. 
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Our proposed mechanism is quantitatively confirmed by 
the values of the Cg hopping integrals from the WF cen- 
tered at Mn(0) site to the one at Mn(l) site (taken from 
off-diagonal matrix element of the Hamiltonian, see Ref. 
[19]): in the optimized AFM-E phase, the hopping inte- 
gral t = 117 meV is enhanced with respect to t = 100 
meV in optimized AFM-A phase (despite the increase in 
the long Mn-0 bond- length after ionic relaxation which 
would counteract the increase in hopping given by the 
increased Mn-O-Mn angle). The hopping in the opti- 
mized AFM-E phase is also increased with respect to t 
= 108 meV in optimized AFM-A atomic coordinates but 
with the AFM-E spin configuration. In parallel, the band 
energy {i.e. on-site energy, taken from diagonal Hamil- 
tonian elements) related to Cg states reduces along the 
AFM-A^ AFM-E adiabatic path with the ionic displace- 
ments, being minimal in the optimized AFM-E phase 
{i.e. 8.98 eV vs 9.03 eV). Therefore, the increase in t 
and decrease in energy in the AFM-E spin-configuration 
is mainly determined by the Mn-O-Mn angle dependence 
and to a less extent by the Mn-0 bond-length. 

We further remark that the movement of the Mn ion 
(cfr Fig. [5|) causes a current whose direction coincides 
with the current by the electron hopping and, there- 
fore, it reinforces the net electronic polarization. It 
should be noted that the atomic displacement (0.04 A) 
is very small compared to the deviation of the Wan- 
nier center (0.26 A). We recall in fact that the atomic 
displacement is just a secondary effect which occurs 
in order to enhance the asymmetric hopping integrals. 
This magnetically-induced mechanism is therefore differ- 
ent from the conventional interpretation of polarization 
in standard ferroelectrics[25j where the atomic displace- 
ment dominates the effect. 

Next, consider the t2g'.dyz orbital which makes a 7r-like 
bonding with surrounding oxygen p orbit als. Because the 
dyz orbital has isotropic symmetry in the plane, the hy- 
bridization with the p orbital depends only on the bond- 
distance. In FiglSb, strong hybridization of dyz orbital 
with 0(2)P-p2 and 0(4)^^-^^ is shown. Similar to the Cg 
orbital case, the electron hops only into Mn(2)-up site 
so that the WFC is displaced in such a way (see green 
arrow in Fig|3)3). Moreover, the atomic displacement 
induced by the Cg orbital - explained above - causes a 
shorter bond length between Mn(0) and 0(3)^^ ions (see 
weight of Mn-t2g WF on 0(3)^p); as a consequence, the 
increased hybridization slightly changes the direction of 
the WFC displacement with respect to the Mn(0)-Mn(2) 
direction. Therefore, in HoMnOa, the anomalous contri- 
butions from Mn-e^ and t2g orbitals almost cancel each 
other along the polar c direction, whereas the 0-p con- 
tribution survives. In order to complete our analysis, we 
therefore show in Fig. [4] the Op WF. What happens is 
qualitatively explained as follows: i) the WFC of 0^:py 
state (forming a 7r-bonding with the same-spin Mn-t2^ 
state) is pulled by the short-distant Mn(l) ion (cfr Fig. 
[4] a); a) the WFC of 0^^:^^ which makes a a-bonding 
with the Mn(0)-e^ state is pulled by the Mn(0) ion in the 



same way (cfr Fig. [4]d); Hi) the WFC of 0^:pz which 
forms a <j-bonding with both Mn(l)- and Mn(0)-e^ states 
is displaced towards the Mn(l) ion but also moves to in- 
crease the Mn-O-Mn angle, with a resulting displacement 
roughly parallel to the c axis (cfr Fig. [4]b); iv) the WFC 
of O^'^'.py which does not form strong bonding with short- 
distant Mn ions shows a small displacement (cfr Fig. [H 
c). 



IV. CHEMICAL TRENDS IN RMnOs 

From the experimental point of view, the magnetic 
trend of i^MnOa was found to be strongly affected by 
the Mn-O-Mn bond angle (j) in the Mn02 plane. From 
i?=La to Gd, where 0=155-146°, the ground state is A- 
type AFM; upon decreasing of 0, there is an intermedi- 
ate "lock-in incommensulate-AF state" which couples to 
a ferroelectric polarization at R=Th and Dy, followed by 
a transition to the E-type AFM observed from i?=Ho to 
Lu, where 0=144-140°. [29] 

When (j) is close enough to 180°, the in-plane FM cou- 
pling in AFM-A phase can be explained on the basis 
of Kanamori-Goodenough rules [28] in the framework 
of Mn-0 "semicovalent bonding" . However, when de- 
creases, the overlap of Mn and O orbitals becomes smaller 
and Kanamori-Goodenough rules are not sufficient to 
explain the complex phase diagram of manganites. In- 
stead, the next-nearest-neighbor antiferromagnetic su- 
perexchange becomes relatively dominant and the E type 
AFM state becomes stable. [29] In the following section, 
we will address quantitatively these general arguments 
and discuss how magnetic interactions are affected by (j). 



A. Structural and magnetic properties 

1. Atomic optimization in the AFM-A phase 

The atomic structure of i^MnOs compounds has 
been investigated by neutron powder diffraction [12] (for 
i?=Pr, Nd, Dy, Tb, Ho, Er, Y), single crystal X-ray 
diffractometry [11] (for i?=Nd, Sm, Eu, Gd) and syn- 
chrotron X-ray powder diffraction measurement [10] (for 
i?=Ho, Er, Tm, Yb, Lu). 

First, in order to discuss the magnetic stability in a 
tiny range of energy, we have used the experimental lat- 
tice parameters, where the volume of the unit cell linearly 
decreases in the series according to the Lanthanide con- 
traction upon increasing the atomic number of R atom; 
we have optimized the internal structural parameters im- 
posing AFM-A spin configuration on Mn-d electrons. 
The reason why we used the AFM-A configuration is 
based on the fact that most of these compounds (where 

= La to Gd) show as ground state the AFM-A spin- 
configuration. 

Second, to discuss ferroelectric properties, the atomic 
structure was optimized with AFM-E spin configuration 



6 



which, as discussed for HoMnOa, breaks the inversion 
symmetry of the system and leads to a magneticahy in- 
duced polarization. The results of this second part will 
be discussed in a later section. 

The optimized structural parameters in AFM-A spin 
configuration, reported in Table Hill and FiglHl give a 
rather regular trend in the series. As for the Mn-0 bond 
length, m and s look rather constant over the series, 
whereas / shows a broad maximum at Gd. Although cal- 
culated inter-planar lengths m show good agreement with 
corresponding experimental values, the in-plane lengths 
/ and s show some deviation (< 10%) with respect to ex- 
periments, so that the JT-distortion is underestimated in 
our calculations. The deviation mainly derives from the 
lack of an exact treatment of the exchange-correlation 
potential (and/or many-body effects) at the LDA/GGA 
level: in our simulations, the rather correlated Mn-d elec- 
trons are described as more delocalized than in real man- 
ganites, so that the JT-distortion is reduced. [30] 

As for the Mn-O-Mn bond angle, it decreases almost 
linearly with the R ionic radius. The inter-plane angle 
becomes larger than in-plane angle for small values of the 
radius or R ion. Moreover, the difference between the 
two Mn-0 distances (/ and s) becomes rather constant 
at = Gd. This implies that the JT-distortion is rather 
"saturated" compared with the GdFeOa-like tilting in a 
Mn-Oe octahedron. 



2. Magnetic Stability 

Using the optimized atomic structure obtained by 
imposing the AFM-A spin configuration, the total en- 
ergy of AFM-A, -C, -G, -E and -E* was calculated 
(FigCJa)). This result shows a very clear trend which 
gives the transition from AFM-A phase to AFM-E phase 
at SmMnOa. Although the optimization imposing AFM- 
A spin ordering enhances the stability of AFM-A, the 
results show that the AFM-E phase is most stable at 
R= Gd, at variance with experimental results report- 
ing that the AFM-A phase is the spin ground-state even 
in GdMnOs. The disagreement with experiments may 
again come from the use of LDA/GGA. Possible improve- 
ments might be obtained upon use of an LDA/GGA+/7 
effective potential [3l| to enhance the localization of Mn-d 
state and reproduce the correct ground state in GdMnOs. 
However, the change of energy is quite sensitive to the 
value of U (which is, by the way, unknown for most of 
the compounds): in our previous calculations [32] , small 
value of U (~ 2.5eV) were found to stabilize AFM-A 
phase with respect to AFM-E in HoMnOs, at variance 
with experiments which show the AFM-E as ground- 
state, [sj, 33] In order to avoid any additional parameters 
and to discuss the trend of the manganites series with- 
out further bias, we don't employ the LDA/GGA+U ap- 
proach and focus on qualitative prediction of chemical 
trends (as a function of the rare-earth ion), which are 
expected to be well reproduced within a bare GGA ap- 




- La Pr Nd Sm Eu Gd Tb Dy Ho Er TmYb Lu- 



1.2 1.15 1.1 1.05 

rR (A) 

FIG. 6: (a) Mn-0 bond distance of of i^MnOs as a function 
of tr (radius of Rare earth ion R^^). Our results are shown 
as round symbols: open and filled circles show long and small 
Mn-0 distances in ac plane, respectively, and the half-filled 
symbols show a middle Mn-0 distance along b axis. Lines are 
plotted for an eye guide, (b) Mn-O-Mn angles of i^MnOs as a 
function of the radius of Rare earth atom. The open symbol 
shows inter-plane Mn-O-Mn angles (with apical O along b 
axis) whereas the closed symbol shows in the ac plane Mn-O- 
Mn angles. In both panels, experimental data are shown for 
comparison: square from [12], diamond from [Tol], left triangle 
from [9.], upper triangle from pj] ]. 



proach. 

Using Heisenberg Hamiltonian with normalized spin 
moment: 

we estimated the super exchange interaction energies Jij . 
From the following six equations, the difference of to- 
tal energy between each AFM phase and the reference 
FM phase is calculated. Then using a least square mean 
method, we obtained four parameters for Jiy. the first- 
nearest-neighbor and second-nearest-neighbor cou- 
pling along a axis in the ac plane, as well as the first- 
and second-nearest-neighbor coupling out of plane J\ 
and J\. With these considered five AFM-configurations, 
the contribution from second-nearest-neighbor coupling 
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TABLE III: Structural parameters optimized with AFM-A configuration in Prima unit cell where the origin is fixed at the 
position on Mn atom. The lattice parameter used for calculations are obtained from Ref.[12] for La, Pr, Nd, Tb, Dy and Ho, 
from Ref.[10] for Er, Tm, Yb and Lu, from Ref.pT,] for Sm, Eu and Gd. 





La 


Pr 


Nd 


Sm 


Gd 


Tb 


Dy 


Ho 


Er 


Tm 


Lu 


a (A) 


5.7473 


5.8129 


5.8317 


5.8620 


5.8660 


5.8384 


5.8337 


5.8354 


5.8223 


5.8085 


5.7868 


b(A) 


7.6929 


7.5856 


7.5546 


7.4770 


7.4310 


7.4025 


7.3778 


7.3606 


7.3357 


7.3175 


7.2959 


c(A) 


5.5367 


5.4491 


5.4170 


5.3620 


5.3180 


5.2931 


5.2785 


5.2572 


5.2395 


5.2277 


5.1972 


4c (x ^ z) 
























X 


0.0508 


0.0654 


0.0714 


0.0785 


0.0829 


0.0836 


0.0846 


0.0856 


0.0864 


0.0867 


0.0870 


z 


0.4907 


0.4866 


0.4849 


0.4827 


0.4816 


0.4808 


0.4806 


0.4805 


0.4798 


0.4799 


0.4800 


Mn 4a (0 0) 


01 4c (x { z) 
























X 


0.4856 


0.4820 


0.4798 


0.4745 


0.4683 


0.4665 


0.4643 


0.4617 


0.4600 


0.4574 


0.4536 


z 


0.5779 


0.5838 


0.5888 


0.5979 


0.6066 


0.6102 


0.6133 


0.6162 


0.6190 


0.6223 


0.6267 


02 8d (x y z) 
























X 


0.3014 


0.3079 


0.3124 


0.3181 


0.3224 


0.3223 


0.3235 


0.3250 


0.3250 


0.3258 


0.3272 


y 


0.0404 


0.0429 


0.0452 


0.0485 


0.0515 


0.0530 


0.0541 


0.0550 


0.0561 


0.0573 


0.0591 


z 


0.2215 


0.2155 


0.2117 


0.2071 


0.2033 


0.2020 


0.2006 


0.1988 


0.1981 


0.1972 


0.1952 



along the c axis is not taken into account and, therefore, 
cannot be determined. 



FM 


: ^ = 4J^" + 2jp + 


2Ji + 


8Ji (2) 


AFM-A 


: E = Uf + 2Jf- 


2Ji - 


8Ji (3) 


AFM-C 


: E = -4Jf + 2Jf 


+ 2Ji 


-8Ji (4) 


AFM-G 


: ^ = -4 J|^" + 2 J|P 


-2J]_ 


+ 8 Ji (5) 


AFM-E 


: E = -2Jf - 2Ji 




(6) 


AFM-E* 


: ^= -2J|P + 2Ji. 




(7) 


As shown in 


FigEfb), consistently 


with 


Kanamori- 



Goodenough rules, [28] the magnitude of ferromagnetic 
R. J|^" represents the sum of two 



J|^" decreases with 



competing interactions (FM coupling due to Cg orbitals 
and AFM coupling due to t2g coupling); the global FM 
behaviour shows that the former dominates. However, 
when (/) decreases, the orbital overlap between Mn-e^ or- 
bitals making a bonding with 0-p orbital is strongly re- 
duced and this, in turn, reduces the FM character of J|^". 
In-plane anti- ferromagnetic J|^"" stays rather constant, 
at variance with what proposed in the model study by 
Kimura et al where J|^"" increases with (j). [29] The inter- 
plane coupling is expected to be AFM, due to its mainly 
t2^-driven character. The fact that it is estimated to be 
weakly FM (< 2 meV) in LaMnOa casts some doubts 
about the validity of a bare GGA treatment for LaMnOa 
(where the AFM-A spin-state is found basically degen- 
erate with the FM spin configuration, obviously at vari- 
ance with experiments). On the other hand, the correct 
spin ground-state is reproduced for most of the i^MnOs 
(with exceptions of i) the above mentioned Gd andii) Tb 
and Dy where we did not attempt to simulate the non- 
collinear spiral arrangements, due to further complexity 
in the simulations). [29] 

Figure [Tl^c) shows the total energy using Eqn.([3]) and 
(|6|): within the mean field approximation, this energy is 



supposed to be proportional to the ordering Neel tem- 
perature Tn- Indeed, the trend is in good agreeement 
with experimental results, showing a steeply decreasing 
Tn with AFM-A in the first half of the series and a rather 
constant Tn with AFM-E in the second half. 

As a summary of this section, an increase of (j) strongly 
reduces e^-derived Jy^ but doesn't change drastically the 
Jj^"" exchange constant, so that the AFM-A phase shows 
a spin transition to AFM-E phase. 



B. The multiferroic AFM-E spin configuration 

1. Structural properties 

In order to discuss the chemical trends of ferroelectric 
properties, we optimized the internal atomic positions in 
the i^MnOs systems by artificially imposing the AFM-E 
spin configuration for all i^MnOs (irrespective of whether 
or not the AFM-E is the spin ground-state). The result- 
ing symmetry (lowered by the spin configuration with 
respect to the AFM-A spin arrangements), allows two 
significantly different values of / and (j) between paral- 
lel and anti-parallel Mn-spins whereas the optimized s 
doesn't change (cfr Fig|8]). This is consistent with our 
mechanism proposed for HoMnOa in the previous sec- 
tion: since the long Mn-0 bond length is mostly affected 
by the ferroelectric AFM-E spin configuration (cfr Fig|8|) 
and it is the one along which the Mn-dz'^ orbital is point- 
ing, it is likely that this latter orbital plays a key role 
in the final ferroelectricity. We also note that, along the 
series, the difference between (jjP and (j)^^ stays rather con- 
stant; since, according to Ref. [s'], this difference is at the 
basis of the polar atomic displacements, this constant be- 
haviour will be relevant in the discussion of polarization 
trends vs R (see below). 
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FIG. 7: (a) Total energy of AFM phases with respect to FM 
phase, (b) Super exchange energy (Jij): The first nearest- 
neighbor in-plane J^^, next nearest-neighbor in-plane J|^^^, 
first nearest-neighbor inter plane J^, next nearest-neighbor 
inter plane J^. (c) Total energy of AFM-A and AFM-E phase 
within mean- field approximation (i.e. sum of Jij)- 



2. Electronic States and Orbital Ordering 

In order to better understand the trends of the ferro- 
electric properties as a function of i?, reported in the fol- 
lowing section, let us briefly discuss the electronic struc- 
ture of the compounds in the AFM-E spin configuration. 

As far as the magnetic moment in the Mn muffin-tin 
sphere is concerned, it is basically constant and equal 
to ~ 3.31 jiiB for all R. However, some differences arise 
when looking at the partial density of states (pDOS) for 
Mn-3(i state along the series (projected in the "local" oc- 
tahedron frame to highlight orbital-ordering), as reported 
in Figure m It is clearly evident that the d^^ and dx2_y2 
states are fully orbitally-polarized. More quantitatively, 
the coefficient of the orbital, obtained by diagonaliza- 
tion of density matrix for five d orbitals, is 0.98 at i?=La 



2.3 
< 2.2 

I ^-^ 

O 

i ' 

1.9 
^ 155 




La 



Pr Nd 



Sm Eu Gd Tb Dy Ho Er Tm Lu: 



1.2 



1.15 „ 1.1 

>-r(A) 



1.05 



FIG. 8: The Mn-0 bond length (upper panel) and Mn-0- 
Mn bond angle (lower panel) as a function of rn in i^MnOs 
structure optimized in AFM-E phase. Here, the superscript 
p (ap) denotes parallel (antiparallel) spin of Mn ions. 
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FIG. 9: The orbital-decomposed partial DOS of Mn-Sd com- 
ponents for i^MnOs : = a) La, b) Nd, c) Ho, d) Lu. Arrows 
indicate: Energy gap (Eg), crystal field splitting (^cf), eg 
bandwidth (w) and on-site exchange interaction energy (J). 
The crystal structure was optimized imposing the AFM-E 
spin configuration. The local frame in the MnOe octahedron 
is chosen as x:middle(inter), y:short, z:long axis. 
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TABLE IV: Hopping integral between Mn-e^ WFs with same 
spin state in optimized AFM-E i^MnOs (meV). 







Nd 


Sm 


Ho 


Lu 


t 


128 


128 


124 


117 


112 




- La 



Pr Nd 



Sm GdTb DyHo ErTm Lu - 



155 



150 



145 



140 



(j) (deg) 



FIG. 10: Ferroelectric polarization calculated by Berry phase 
method (BP), Wannier function method (WF) and point 
charge model (PCM) as a function of in AFM-E i^MnOs 
(see text for details). The in-plane Mn-O-Mn bond angle (f) is 
averaged between 0^ and (/)ap. 



and 0.96 at R=Lu. This implies that the orbital ordering 
is already "saturated" at LaMnOs. 

A progressive distortion of the structure is shown to 
increase the energy gap Eg and to decrease the width 
of the e^-band {w). Within a tight-binding framework, 
w is proportional to the hopping integral t. Indeed, we 
estimated the hopping integral for selected compounds 
along the series and found that t increases with the ionic 
radius of the R ion (see Table ITVj) . 

Here, we would like to point out that the DOS doesn't 
change significantly before and after FE atomic displace- 
ments, at variance with standard FE (such as BaTiOs), 
where FE atomic displacements are accompanied by a re- 
hybridization of filled 0-p and empty cation d states. Q 
Here, we recall that the Mn-d state is well localized and 
the character of the Mn-0 bond is rather ionic in nature 
so it is not expected to undergo drastic changes upon de- 
velopment of polarization as far as the bonding properties 
are concerned. 



3. Ferroelectricity and link to lattice, spin and orbital 
degrees of freedom 

Similarly to the case of HoMnOs, we have calculated 
the polarization according to the PCM, BP and WF ap- 
proaches. As shown in FigjlOl the polarization calculated 
by BP and WF are consistent within less than 10 %: this 
difference (in addition to possible numerical uncertainties 
deriving from the use of a different basis set and potential 



TABLE V: Polarization (/iC/cm^) induced by the displace- 
ment of Wannier function center (WFC) from the ionic po- 
sition in AFM-E i^MnOs. Only up-spin contribution is re- 
ported for each WF: Mn-e^, Mn-t2^ and 0-p. In the two final 
lines, we show the sum of the contributions (for spin-up only) 
and the total (for spin-up and spin-down). 



p 


= La 


Nd 


Sm 


Ho 


Lu 


Mn-Cg 


-4.76 


-3.48 


-3.01 


-2.67 


-2.61 


Mn-t2g 


1.36 


1.73 


2.02 


2.42 


2.47 


0-p 


-0.96 


-1.00 


-1.17 


-1.51 


-1.55 


total (up) 


-4.36 


-2.75 


-2.16 


-1.76 


-1.69 


total (up+down) 


-8.72 


-5.50 


-4.32 


-3.52 


-3.38 



treatments in VASP and FLEUR codes for BP and WF 
approaches, respectively) comes from 0-2s, R-bs, R-bp 
contributions and shows that these latter contributions, 
as expected, are small compared to Mn-d and 0-p con- 
tributions. 

What is truly remarkable is that the polarization (both 
from BP and WF) shows a rapid increase upon decreas- 
ing of (j) (Pbp approaches 12 (iC /cm? in an hypothetical 
AFM-E LaMnOa), whereas PpcM is to a large extent con- 
stant. This latter trend is consistent with the constant 
behavior of 0^ — (j)^^ shown in Fig|8] and with what pre- 
viously discussed in Refs. [3, 7]. Moreover, the difference 
between Pbp /Pwf and PpcM (which represents a purely 
electronic contribution) is much bigger than the contri- 
bution coming from atomic displacementsin the entire 
series. This confirms the different nature of the magnetic 
origin of ferroelectricity in manganites with respect to 
standard "proper" ferroelectrics. 

Our results further suggest that orbital ordering is nec- 
essary for the rising of the polarization; however, being 
constant along the series, is not responsible for the trend 
of FE polarization, which, on the other hand, shows dra- 
matic changes as a function of R. Consistently with what 
previously discussed for HoMnOs where t is suggested to 
play a important role, the dramatic increase in the asym- 
metric Cg hopping is responsible for enhancing P in less 
distorted manganites (cfr Tab. [Vl where we report the 



different contribution of Mn-e^, Mn-t2^ and 0-p to the 
total Pwf)- We show that, for every i?, the Mn-e^ and 
0-p contributions show an opposite sign with respect to 
the Mn-t2^ term with regular trends along the series. In- 
deed, remembering that the hopping integral t depends 
on the bond length and bond angle, one expects smaller 
/ and larger (j) to enhance t between Cg orbit als, whereas 
smaller s to enhance t between t2g orbitals; this is con- 
firmed by looking carefully at tablelVland Fig 15] and their 
implications for the values of P (cfr. Fig. [K 



V. CONCLUSIONS 

Ferroelectricity, recently proposed for the collinear 
AFM-E magnetic phase in orthorhombic HoMnOa, is 
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here explained microscopically from first-principles via 
a careful Wannier function analysis. We show that the 
asymmetric electron hopping of orbitally-polarized Mn- 
Cg states is the key ingredient for the rising of polariza- 
tion. At variance with proper ferroelectrics, in HoMnOa 
the purely electronic contribution (due to Wannier func- 
tion centers which are largely displaced with respect to 
ions) dominates the polarization, with respect to the con- 
tribution to P coming from ionic movements. However, 
the net polarization along the c axis is the result of a 
delicate balance of different contributions, such as, for 
example, the opposite signs of the Mn Cg and t2g contri- 
butions to P. 

In addition, extensive ab-initio calculations have been 
performed for the i^MnOa systems, focusing on the link 
between ferroelectricity and the spin, orbital and lattice 
degrees of freedom in the aim of identifying chemical 
trends along the series. In summary, the results are sum- 
marized as: i) ferromagnetic J^^ increases with the Mn- 
0-Mn angle (/), whereas antiferromagnetic J^^^ is rather 
constant with 0; ii) P^^ increases with (j) and decreases 
with /, whereas Pt^g decreases with s and partially can- 
cels Pgg • ^^^) orbital ordering, a needed ingredient in the 
rising of P, is however saturated along the series and does 



not influence the trend of polarization as a function of R. 

Therefore, the main message is that one can expect 
high ferroelectric polarization with large (j) and small / 
in the AFM-E phase; however, such large (j) would si- 
multaneously increase the ferromagnetic J^^^ therefore 
stabilizing the centrosymmetric AFM-A spin configura- 
tion instead of the polar AFM-E phase. The problem is 
then how to increase P, still keeping the AFM-E phase as 
stable spin-state {i.e. via strain, pressure, alloying, etc): 
we hope that our findings will be helpful to answer this 
question. 
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